The role of screening in the density functional applied on transition metal defects in semiconductors 
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We study selected transition metal related point defects in silicon and silicon carbide semiconductors by a 
range separated hybrid density functional (HSE06). We find that HSE06 does not fulfill the generalized Koop- 
mans' Theorem for every defect which is due to the self-interaction error in the functional in such cases. Restor- 
ing the so-called generalized Koopmans' Condition with a simple correction in the functional can eliminate this 
error, and brings the calculated charge transition levels remarkably close to the experimental data as well as to 
the calculated quasi-particle levels from many-body perturbation theory. 
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Tractable solutions for the problem of many-electron sys- 
tems are highly needed in materials science and under intense 
research. One key example of such problems is the treatment 
of small imperfections in single crystals, i.e. point defects in 
semiconductors. Point defects can completely change locally 
the chemical bonds of the host crystal, thus ab initio methods 
are needed in order to determine their geometry, electronic 
structure, ionization energies or optical excitations. The most 
challenging type of point defects is the transition metal (TM) 
related defects in traditional semiconductors where d-orbitals 
tightly localized on the TM atoms and sp 3 hybrid orbitals are 
both present in the system. It is known from , e.g. studies of 
Mott-insulators that d-electron systems may be highly corre- 
lated whereas the electronic structure of traditional semicon- 
ductors can phenomenologically be well-described by inde- 
pendent particle theories. Simulations of TM defects in semi- 
conductors are therefore a complex problem, thus it is an obvi- 
ous system to investigate the predictive power of the ab initio 
methods. 

Density functional theory (DFT) is the most widespread 
technique for first principles calculation in condensed matter 
physics. Particularly, DFT has been proven to be extremely 
powerful tool to study defects in semiconductors [1]. The 
success of the DFT calculations is based on well-developed 
approximate (semi)local functionals [2-5 ] that made it pos- 
sible to study relatively large systems at moderate computa- 
tional cost with a surprisingly good accuracy [1]. The suc- 
cess of the commonly used (semi)local functionals might be 
unexpected as they suffer from the self-interaction error Q 
which results in the underestimation of the band gap of semi- 
conductors. The predictive power of (semi)local functionals 
is thus restricted for ionization energies of defects in semi- 
conductors, and can even fail to describe the nature of the 
host semiconductor, or the defect state in the semiconductor 
correctly in pathological cases (SHU- By mixing a non-local 
Fock exchange into the density functional, i.e., by using hy- 
brid density functionals, one can restore the band gap of host 
semiconductors via tuning the mixing parameter l9i{T3ll. Al- 
ternatively, by introducing a range separated hybrid density 
functional (HSE06) it is found that with a fixed mixing pa- 
rameter and the range of separation, the band gap of many 



semiconductors with sp 3 hybrid orbitals can be well repro- 
duced lfT4l . and in addition, the structural parameters, redox 
reaction energies, and formation energies of transition metal 
compounds can also be well accounted by HSE06 [15]. This 
brings a hope that HSE06 is able to predict the ground state 
lfT6l and ionization energies for TM related defects in semi- 
conductors which contain both d and sp 3 orbitals. However, a 
consistency of the approach requires the use of the same mix- 
ing and range separation parameters for the host and impurity 
states. It has been observed that this requirement makes it 
difficult to describe band structure and defect states simulta- 
neously 071(151 . 

In this Letter we study TM defects in silicon (Si) and AH 
silicon carbide (SiC) semiconductors by HSE06 functional 
using plane wave large- supercell calculations, for those ex- 
perimental data about their structure and ionization energies 
are available. We show that HSE06 can qualitatively fail in 
such complex systems due to insufficient screening of the 
Coulomb interaction between localized d-electrons, leading 
to incomplete cancelation of the self-interaction error which 
manifests as disobeying the generalized Koopmans' Theorem 
(gKT) O [Tt] H21 EQl - We suggest a correction scheme for the 
hybrid functional which fulfills the conditions of gKT. The 
corrected functional allows for a simultaneous description of 
localized and extended states, is first-principles in nature, and 
brings our theoretical results remarkably close to the exper- 
imental data as well as to the results obtained by means of 
many -body perturbation theory. Our work underlines the im- 
portance of testing a fulfillment of gKT for the applied func- 
tionals on a given many-electron system. 

We selected interstitial Fe in Si (FeO, and titanium (Tisi), 
vanadium (Vsi) and chromium (Crsi) substituting Si- site in 
AH SiC. For all of these defects, the structure and the cor- 
responding ionization energies are well established from ex- 
periments ll2TU23l [25l . In addition, we investigate tungsten 
(W) in AH SiC because W-related ionization energies have 
been detected recently [24]. In this case, the origin of the 
W-related centers are not ambiguously identified. According 
to our previous study l26lL both Si- substitution site (Wsi) and 
asymmetric split vacancy configuration (Wasv) may be stable, 
thus both defects are considered. For the electronic structure 
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FIG. 1 . (Color online) Measured and theoretically predicted charge 
transition levels of transition metal impurities in AH SiC and Si 
hosts. In all cases the experimental results [21-24] are represented 
with dark gray (first) columns where the thick horizontal line repre- 
sents the charge transition levels. The result of HSE06 calculations 
and our correction method with satisfying gKT are represented with 
red (second) columns and light green (third) columns, respectively, 
where thin horizontal lines provide the calculated charge transition 
levels. The bars indicated by dashed lines around the calculated 
charge transition levels show the inherent uncertainty of the calcu- 
lation. Chromium needs special consideration regarding the (0|-) 
transition level, see text. 

calculations we use VASP package l27l with plane wave basis 
(cut-off: 420 eV) and projector augmented- wave (PAW) l28l 
method for the ion-cores. In the case of metal atoms we uti- 
lize small core PAW projectors. In order to model AH SiC and 
Si host we use a 576-atom and 512-atom supercells, respec- 
tively. In the large supercells the Brillouin zone is sampled at 
T-point. In the case of charged defects the size dependence of 
total energies is eliminated by charge correction l29l with the 
order of magnitude 0.1 eV in our supercells. The geometry of 
the defects is optimized unless the forces acting on the atoms 
are less than 0.01 eV/A. 

In order to examine the performance of the HSE06 screened 
non-local hybrid density functional on the selected transition 
metal defects in semiconductors we calculate the adiabatic 
ionization energies or charge transition levels of these defects, 

e (q\q') = (E q , + £ corr (</)) - (E q + (q)) - s CBM (1) 

where E q is the total energy and E con (q) is the charge correc- 
tion of charge state q and Eqbm is the conduction band min- 
imum (see Ref. 1 for more details). The HSE06 functional 
gives the proper description of the band structure of the semi- 
conductor host. The band gap error is within 0.1 eV for the 
pristine AH SiC and Si host which causes uncertainty in the 
calculated e (q\q f ) of about 0.1 eV. By comparing the experi- 
mental data with the results of HSE06 calculations it is appar- 
ent from Fig.[T]and Table[T]that HSE06 cannot predict the cor- 
rect charge transition levels for all the defects. The magnitude 
of the discrepancies depend on the TM as well as the charge 



TABLE I. (Color online) Experimental values of charge transition 
levels (e e xp), deviations of the theoretically predicted values from the 
experiment for transitional metal related point defects, Tisi, Vsi, Crsi 
and Wsi in AH SiC and Fei in Si. In the calculations HSE06 func- 
tional and a corrected HSE06 functional (see text) are used. The 
finite values of £?nk clearly indicates the presence of self-interaction 
error in the HSE06 functional (A£hseo6). When small value of £^nk 
is obtained (.Enk) with the help of additional correction functional 
(V w ) with a strength parameter w, the discrepancy (A£hseo6+v w ) goes 
below 0. 1 eV which is the uncertainty in our calculation. Light red 
(dark green) numbers represent the case when the calculated charge 
transition level is out of (within) the error bar. The unit of the data is 
eV in all columns. 
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a Ref. |22], b Ref. (24), c Ref. (21] 



transition level (TableJI]). Particularly, the calculated ionization 
energies of Cr strikingly differs from the experimental values, 
and are qualitatively wrong. 

Let us demonstrate that the failure of HSE06 for some de- 
fects originates from the incorrect treatment of atomic-like 
d-orbitals of TM atoms. It is well-known that self-Hartree 
and the exchange potentials do not cancel each in (semi)local 
DFT functionals. In (semi)local functionals the resultant self- 
repulsive potential lower the localization degree of the states 
and cause the spurious occupation dependence of the Kohn- 
Sham energies as well as the spurious, generally convex, cur- 
vature of the total energy with respect to fractional occupa- 
tions 0[l9][2O]l. The opposite behavior is characteristic for 
Hartree-Fock (HF) method. By mixing HF and semi-local ex- 
change, the self interaction error is reduced in hybrid func- 
tionals. Unfortunately, it is not eliminated completely. 

Previous studies show ||7l H21 H21 EQl that one of the quanti- 
tative manifestation of the self-interaction error in DFT func- 
tionals is the discrepancy between the Kohn-Sham eigenvalue 
of the highest occupied state and the corresponding ioniza- 
tion energy. This energy difference is usually called as Non- 
Koopmans' energy, 

E^k = €n - Hi = e N — (E N — En-i) (2) 

where Sn is the Kohn-Sham energy of a localized state in 
N electrons system, Ej is the ionization energy of the sys- 
tem which equal to the difference of the total energies of the 
TV electron, En, and TV — 1 electron system, E^-i- In HF 
theory the Koopmans' Theorem states that the single particle 
energies equal to the ionization energy in case of every oc- 
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cupied states, while in DFT this statement is only valid for 
the highest occupied state with an exact exchange-correlation 
functional. In the case of exact functional £^nk = condition 
should be fulfilled, the so-called generalized Koopmans' Con- 
dition (gKC). Due to the Janak's and Slater's theorems, other 
important features of the exact functional is simultaneously 
remedied: linear behavior of the total energy and constant be- 
havior of (highest occupied) single particle level with respect 
to fractional occupation are also fulfilled. The last condition 
explicitly shows that the single particle state does not suffer 
from any occupational dependent potential, thus represents a 
self-interaction free functional. 

In Table U one can see the difference in the calculated 
HSE06 results and experimental data as well as the calculated 
Non-Koopmans' energy for the defect state that gets occu- 
pied through the charge transition. In the calculation of E^k 
we considered the two charge states of the defect with fixing 
the geometry as found in the first charge state. The results 
clearly show that this quantity correlates with the discrepancy 
between HSE06 results and experimental data for the charge 
transition levels. The obvious conclusion of this finding is 
that the errors are originated from the spurious self-interaction 
which still remains in the HSE06 functional. 

In HSE06 functional, the mixing parameter and the range 
of separation play the role of effective screening of Coulomb 
interaction [30-32], and their variation in practice interpo- 
lates between the substantial underestimation of the localiza- 
tion effects in (semi)local functionals and their overestima- 
tion in Hartree-Fock theory. But the strength of the screening 
should itself depend on the localization degree of the states. 
In the case of TM impurities the strong repulsive interaction 
between the electrons stems the extension of the atomic-like 
electron states. In such a strongly correlated system, the lo- 
calized states are more favorable than the extended one. Due 
to the large localization degree of the d- states with respect 
to that of host sp 3 hybrid orbitals, in general it should not 
be possible to describe both subsystems with the same values 
of parameters in a hybrid functional. Lany and Zunger ifTTl 
demonstrated that one should not expect that adjusting param- 
eters of hybrid functionals would automatically give accurate 
results for both, host band gaps and defect levels, and pointed 
out homogeneous screening of the Fock exchange in hybrid 
functionals as serious simplification. Still, it is desirable to 
have an option to describe the orbitals with different degrees 
of localization in calculations. 

At this stage, it is important to mention that, a systematic 
study of the role of nonlocal exchange in the electronic struc- 
ture of correlated oxides by F. Iori et al. l32l demonstrates 
that the default values of parameters of HSE06 functional sys- 
tematically overestimate band gaps in these materials, indicat- 
ing insufficient screening of Coulomb interaction. Negative 
values of £?nk in Table [i] also indicate that the screening of 
Coulomb interaction is insufficient |[T7l . 

In order to prove the self-interaction error of HSE06 func- 
tional, we carry out sophisticated parameter free GWo cal- 
culations 1 33] on defective supercells with 128 atoms and 



TABLE II. The position of the highest occupied single particle levels 
of transitional metal defects in the gap with respect to the conduction 
band edge (CBM) in eV unit. The table shows the results of calcu- 
lations with GWo method, HSE06 functional and corrected HSE06 
functional (HSE06+V W ). In the case of Crsi(0) the corresponding 
single particle state falls in the valence band shown as "VB", thus 
GWo started far from the true ground state of the system. When Non- 
Koopmans' energy (£?nk) is not negligible the correction method as 
well as the GWo method shift the levels up with approximately the 
same value. This result confirms the presence of self-interaction er- 
ror in the HSE06 functional and validate our corrected functional. 
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HSE06 relaxed geometry. We apply 1344 bands in the cal- 
culation of the response function, and the Brillouin zone is 
sampled with 2x2x2 Monkhorst-Pack k -point set (34). The 
starting wave functions are obtained from HSE06 calculation, 
then the Green-function G and the wave functions are self- 
consistently updated while keeping the screened Coulumb- 
interaction W fixed. We find that four iterations are suffi- 
cient to reach the self-consistent quasi-particle levels within 
0.05 eV. The observed quasiparticle correction (see Table |n] ) 
confirmed the over-localization error of the HSE06. 

Considering £? NK is an appropriate measure to construct 
a self-consistent correction method to eliminate the self- 
interaction error in the functional (71 [T9lL we propose a 
scheme that does not use any empirical parameters to coun- 
teract the error due to insufficient screening of the Coulomb 
interaction between d-electrons localized at TM related de- 
fects, which is apparently present in HSE06 functional. 

In our correction technique the occupation dependent po- 
tential 

V£ TO = |(l-2r4) (3) 

is applied on the d-orbitals together with the screened non- 
local HSE06 functional. In Eq. ^ w is a parameter for the 
strength of the potential, which physical meaning will become 
clear below, and which is determined self-consistently by sat- 
isfying gKC. We emphasize that the resultant functional is ab 
initio in the sense that it is not fit to any empirical parameter. 
We note that the form of Eq. ^ is equivalent to Dudarev's 
implementation of LDA+U method |35], and therefore it is 
straight-forward to use the suggested scheme with existing 
first-principles package like VASP. The occupation number 
in VASP implementation is determined as projection of 
the wave function of the system on the spherical harmonics 
Im on site of atom /, where Z = 2 in our case. 



3 



The results are summarized in Table [I] One sees that gKC 
could be fulfilled with great accuracy (E^ K ) in all cases ex- 
cept for Fei in Si and Crsi : (0|-) in SiC. In case of iron E^k 
has an extrema and does not reach zero. We observed that 
the particular defect state is not so well localized on the Fe 
d-orbital in Fei defect, thus our correction scheme is not suf- 
ficiently effective. Similar effect can be found in case of neg- 
atively charged Crsi where in the ground state three electron 
occupy the e and a\ levels with parallel spins. The highest 
occupied a\ level is also less atomic-like compared to the e 
level. In Table [i] the Crsi:(0|-) charge transition level is cal- 
culated as sqw + £; relax + E cox \ where eqw is the highest 
occupied quasi-particle level in the negative charge state, and 
E TQlax is the relaxation energy of the defect due to ionization 
(36). Note that achieving gKC with a certain accuracy brings 
the corrected single particle level to be independent from the 
occupation number with the same accuracy indicating the self 
interaction free description of the defect state. The charge 
transition level was obtained after the geometry relaxation of 
system with the use of the corrected HSE06+V W functional 
shown in the last column of Table [l| Apparently, this correc- 
tion method can bring the theoretical values closer to the ex- 
perimental values in every case. When gKC is accurately sat- 
isfied the calculated levels approach the experimental levels 
within 0.1 eV which is the expected uncertainty in our cal- 
culation. It also can be seen that when gKC is fulfilled the 
position of the defect state in the gap is close to the corre- 
sponding quasi-particle level calculated by the GWo approx- 
imation of the many-body perturbation theory. We note here 
that the formation energies of the defects are also affected by 
our correction scheme which is discussed in the Supplemental 
Material l36l . 

Considering the strength of the required correction poten- 
tials an important remark can be made. In all the cases when 
the correction is needed, the parameter w is negative repre- 
senting a repulsive potential for occupied states. This ob- 
servation confirms that HSE06 functional over-localizes the 
defect state of the TM impurities, while this spurious effect 
is eliminated with our technique. In more details this is il- 
lustrated in Fig. [2jb) where the difference of the spin den- 
sity of HSE06 and HSE06+V W calculation is shown. In the 
case of TM impurities, the HSE06 functional clearly repre- 
sents an over-correction with respect to standard (semi)local 
DFT functionals. 

In order to understand the physical meaning of parameter w 
in Eq. ([5]), we will make use of the similarity between the way 
how hybrid functionals and LDA+U approach (35j [37] |38) 
correct the self-interaction error, pointed out by several groups 
l32l . For the case of correlated electrons it becomes par- 
ticularly apparent within the so-called exact exchange for 
correlated electrons method by Novak and co-workers l39l . 
LDA+U functional was introduced in order to take into ac- 
count the effects of the strong correlation which is usually not 
properly described by (semi)local functionals l35ll37l[38l . 

The positive value of the Hubbard parameter U represent- 
ing the strength of the screened Coulomb potential makes the 
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FIG. 2. (Color online) (a) The calculated electron density (with an 
isosurface value of 0.05) of the unpaired d-orbital of the neutral vana- 
dium substitutional defect in 4H SiC. (b) 2D plot of the differences 
between two calculated charge densities using HSE06 functional and 
the corrected HSE06 functional around the vanadium atom. The pos- 
itive finite density clearly represents the overlocalization error of the 
HSE06 functional. 



atomic like states more favorable and increase the localiza- 
tion of the states on atomic orbitals. In practice the screening 
of U is achieved by a reduction of the large bare Coulomb 
value, and the actual strength of the on-site electron repulsion 
can be either determined empirically l40ll4TlL or evaluated ab 
initio (42l [43). As a matter of fact, in our case of the TM 
defects the Coulomb interaction between localized electrons 
(Fig. [2]) appears to be insufficiently screened in HSE06 func- 
tional due to its homogeneous nature and parameters adjusted 
to describe more extended states. In LDA+U formalism this 
error could be corrected by a reduction of U parameter. Em- 
ploying the above mentioned similarity between the hybrid 
functionals and LDA+U method, we see that the action of oc- 
cupation dependent potential in Eq. ^ on localized orbitals at 
impurity sites effectively increases the screening of Coulomb 
interaction between electrons occupying these orbitals. 

Finally, we show the predictive power of our methodol- 
ogy demonstrated on W-related defects in AH SiC. We find 
that Wsi defect has two charge transition levels and these are 
very close to the experimental values. In the case of Wasv 
defect HSE06 again suffers from the self-interaction error, 
^nk=-0.18 eV, and correction is needed in the functional with 
w=-2A eV. It is worthy to notice that this value differs from 
the value needed for Wsi defect which shows that the self- 
interaction error for TM defects may be different in various 
defect configurations. With our correction method the (0|-) 
charge transition level of Wasv defect is at £cbm-0.64 eV 
where no W-related charge transition level was found in the 
experiments (see Table [I]). This indicates that Wsi is associ- 
ated with the measured W-related center l24l . thus we show 
the calculated ionization energies of Wsi in Table |T| 

In summary, our theoretical investigation revealed that the 
HSE06 hybrid functional may overlocalize the defect states 
arisen from transition metal defects. The reason of this er- 
ror is the self-interaction that is increased due to the insuffi- 
cient screening of Coulomb interaction between electrons lo- 
calized at the transitional metal defects. By invoking the gen- 
eralized Koopmans' Condition we were able to indicate the 
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presence of self-interaction error in the calculation. Further- 
more, we suggest a correction method which can eliminate 
the self-interaction error of hybrid functionals in the case of 
atomic-like defect states. With this technique all the calcu- 
lated charge transition levels become accurate within 0.1 eV. 
In addition, the corrected Kohn-Sham levels are close to those 
obtained by the sophisticated but computationally demanding 
GW -method. 
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1. Calculation methods 

For the electronic structure calculations we use VASP package [Rl] with plane wave basis (cut-off: 
420eV) and projector augmented- wave (PAW) [R2] method for the ion-cores. In the case of metal 
atoms we utilize small core PAW projectors. In order to model 4H SiC and Si host we use a 576-atom 
and 512-atom supercells, respectively. In the large supercells the Brillouin zone is sampled at T-point. 
In the case of charged defects the size dependence of total energies is eliminated by charge correction 
[R3] with the order of magnitude 0.1 ~eV in our supercells. The geometry of the defects is optimized 
unless the forces acting on the atoms are less than O.OleV/A. 

GW calculations are carried out in 128-atom supercell where HSE06 relaxed geometry is applied in 
defective supercells. We apply 1344 bands in the calculation of the response function, and the Brillouin 
zone is sampled with 2x2x2 Monkhorst-Pack £-point set [R4]. The starting wave functions are obtained 
from HSE06 calculation, then the Green-function G and the wave functions are self-consistently 
updated while keeping the screened Coulumb-interaction W fixed. We find that four iterations are 
sufficient to reach the self-consistent quasi-particle levels within 0.05 eV. 



2. Relative stability and charge transition level diagram 

Our method is designed to correct a localized atomic like single particle orbital by recovering its real 
charge density. Accurate charge transition levels can be obtained from the sum of the vertical ionization 
energy obtained from the generalized Kohn-Sham level and the geometry relaxation energy due to 
ionization (See Fig. SI and Table SI in paper). The necessary correction on the Kohn-Sham level 
depends on the given transitional metal (TM) as well as on the total occupation number of the d-like 
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states of this TM-atom, or in other words, on the charge state of the defect. Due to the varying 
additional screening functional the total energy will be not "unique". The main difficulty may arise 
from this method when the total energies of different configurations of the given TM impurity are 
compared in order to determine their relative stability. 

To overcome this difficulty, it is possible to choose a common screening parameter w c for two or 
more neutrally charged defects. A rational choice is the average of the parameters w obtained for 
different configurations of the given TM impurity, however it is unknown how sensitive is the relative 
stability to the choice of the common parameter w c . 

We studied the formation energies of the experimentally most relevant TM defects in 4H-SiC. We 
considered the substitutional defect (TM s at h site and the asymmetric split vacancy defect (ASV), 
consisting of TM S i and a carbon vacancy (TM S i-Vc ) at hh site when TM is V or W. 
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(a) (b) 

FIG. SI. (Color online) The total energy (E to t) of substitutional and ASV 
defects of neutral (a) vanadium (V) and (b) tungsten (W) is plotted as a 
function of the applied screening potential (w). 

The total energies of the two types of TM-defect with respect to screening parameter w can be seen 
in Figs. SI (a) and (b). Total energies depend linearly on the parameter w and, in addition, the slope of 
the curves is close to each other, for example 0.912 for V-ASV and 0.885 for Vsi. Due to this behavior 
the relative stability is approximately (within -0.1 eV) independent on the choice of the common 
screening parameter w c and is close to results of HSE06 calculation. In critical cases (e.g. when the 
difference in total energies is close to 0.1 eV), the hierarchy may change as a function of w parameter. 
W impurity is such an example in 4H SiC. 

The obtained screening parameter and the corrected charge transition levels for the considered 
defects are summarized in Table SI. 
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TABLE SI. The optimized screening parameter w that fulfill the gKC and the 
corresponding corrected charge transition levels with respect to the conduction band 
edge (AEhseo6+Vw) of substitutional and ASV defects of V and W. 



defect 


w 


AEhsE06+Vw 


[eV] 


[eV] 


V Si : (+10) 


-2.7 


- 1.69 


V Si : (01-) 


-2.2 


-0.94 


V-ASV : (+10) 


-4.6 


- 1.47 


V-ASV : (01-) 


-4.2 


- 1.29 


W si : (+10) 


-1.2 


- 1.43 


W si : (01-) 


0.0 


-0.12 


W-ASV : (+10) 


-2.6 


-2.44 


W-ASV : (01-) 


-2.4 


-0.64 



The chosen common parameters and the corresponding relative stability at Si-rich condition can be 
found in Table SII. With this data it is possible to draw the transitional level diagram which shows the 
relative stability with respect to the Fermi-level (Figure S2). We note that the HSE06+V W functional 
with the use of common w c parameter produces superior results in comparison to conventional 
calculations using HSE06 functional. 



TABLE SII. The chosen common screening parameters w c for the pair of substitutional 
and ASV defects of V and W. The relative stability (E m \) of neutral defects is shown as 
obtained by HSE06 and HSE06+V W functional. 





W c 




£rel,HSE06+V 


defect pairs 


[eV] 


[eV] 






[eV] 


V Si - V-ASV 


-3.65 


1.82 


1.87 


W Si - W-ASV 


-1.9 


0.03 


-0.15 
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FIG. S2. (Color online) Charge transition level diagram of two different 
types of defect of vanadium (V) (a) and tungsten (W) (b) at Si-rich condition 
as obtained by HSE06 and HSE06+V W functionals (see text). Black vertical 
lines represent the position of the experimental charge transition levels. The 
formation energy of the substitutional defect was set to zero as a reference 
energy in both cases. 



3. Partial density of states 

Partial density of states of neutrally charged substitutional vanadium in 4H-SiC is plotted in Fig. S3. 
Calculations are carried out using both HSE06 and HSE06+V W functionals. In the case of vanadium 
(V 4+ ) there is only one electron to occupy the lowest double degenerate d-like state (e) that splits into 
two non-degenerate levels in Cih symmetry. The major spin channel is chosen to be the spin-up 
channel. In the major spin channel both states appear in the gap, the lowest is occupied and the other 
one is unoccupied as the Fermi-level is at 0.37 eV and 1.27 eV in HSE06 and HSE06+V W calculations, 
respectively. The states are slightly hybridized with the s and p orbitals of the host semiconductor. The 
other d-states, however, are so strongly hybridized with the s and p orbitals that they lose their atomic 
character. This is particularly apparent for the occupied states below the Fermi-level [see Fig. S3(c)] 
while empty resonant d-states may occur in the conduction band. The effect of the correction potential 
on those states is almost negligible but very strong on the states occurring in the band gap. 
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FIG. S3. Partial density of states substitutional vanadium in 4H-SiC 
calculated by (a) HSE06 functional and (b) and HSE06+V W functional close 
to the Fermi-level; (c) by HSE06 and HSE06+V W functional in wide energy 
range. The energy of the valence band edge was chosen to zero as reference. 
Up(down) labels the spin-up(down) channel. The spin-up channel is the 
majority spin channel. The calculated Fermi-level is at 0.37 eV and 1.27 eV 
in HSE06 and HSE06+V W calculations, respectively. 



4. Calculation of the electronic structure of chromium 

We mention here that chromium (Cr) defect needs special consideration in determining w. In the ground 
state of the negatively charged Crsi defect a double degenerate e level and an a\ level appear in the gap, 
originated from the split d-orbitals of Cr in the crystal field of 4H-SiC with Cj v symmetry. In the neutral 
charge state two electrons occupy the e state while the a\ state is empty, forming 5=1 spin state. Due 
to the strong correlation of the d-orbitals S - 3/2 state is formed in the negatively charged state where 
three electrons occupy the e and a\ levels with parallel spins. In the negatively charged state, our 
correction method acts on the highest occupied a\ level which is less atomic like compared to the e 
level. Our correction method makes the atomic-like defect states more or less favorable therefore 
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changes the charge density in a way that changes the character of the a\ level. In Table I of the main text 
the Crsii(OI-) charge transition level is calculated as 6bw + £ relax +£ corr , where 6bw is the highest 
occupied quasi-particle level in the negative charge state, and £ relax is the relaxation energy of the defect 
due to ionization. In the case of the Cr S i:(-l2- ) transition, the a\ level becomes empty and the e level 
gets fully occupied by ionizing from (-) to (2-) states. Here, we set S = 1/2 instead of the energetically 
favorable S = 3/2, for the initial negative charge state in order to be able to study of the effect of the 
occupation of e level and to correct its level via satisfying gKC with our correction functional. 



References 

[Rl] G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 (1996) 
[R2] P. E. Blochl, Phys. Rev. B 50, 17953 (1994) 

[R3] C. Freysoldt, J. Neugebauer, and C. G. Van de Walle, Phys. Rev. Lett. 102, 016402 (2009) 
[R4] H. J. Monkhorst and J. K. Pack, Phys. Rev. B 13, 5188 (1976) 



6 



SiC CBM 
3.26 eV 



Charge states 

Columns: 

1. ■ Exp. 

2. ■ HSE06 

3. □ Condition of gKT 



Si CBM 
1.17 eV 



VBM 
eV 
























+ 































+ 



+ 



+ 



+ 



Fei Ti Si V Si 



Cr Si 



W Si 



